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Abstract 


By including the contribution of the thermal cloud to the Lagrangian of the 
condensate of a Bose gas, we extend the time-dependent variational method 
at zero temperature to study temperature-dependent low collective excitation 
modes. A Gaussian trial wave function of the condensate and a static distri- 
bution density of the thermal cloud are used, and analytical expressions for 
temperature-dependent excitation frequencies obtained. Theoretical results 


are compared with measurements in the JILA and MIT experiments. 
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With the development of techniques to trap and cool atoms, Bose-Einstein condensation 
(BEC) has been observed directly in dilute atomic vapors [1-8]. The new experimental 
achievements have stimulated great interest in the theoretical study of inhomogeneous Bose 
gases. The collective excitations of trapped Bose gases, which describe their dynamics and 
transport properties, have been a focus of recent experimental and theoretical studies. Low- 
lying collective excitations over a range of temperatures ИН and higher-lying modes В] 
have been measured in recent experiments. 

In the regime with no detectable non-condensate fraction, theoretical predictions based 
on a mean-field description of weakly interacting bosons at zero temperature are in excellent 
agreement with experimental data [HIT]. At zero temperature, a Bose-Einstein conden- 
sate is described by the nonlinear Schródinger equation (NLSE), i.e. the Gross-Pitaevskii 


equation |12], 


DU. Ла Anh? a 


ih РФ, (1) 


m 


where m is the atomic mass, U(r) the trapping potential, and a the S-wave scattering length. 
Characterizing the evolution of the macroscopic wave function ф of the condensate, this 
nonlinear equation provides a fundamental means for studying the ground state, dynamical 
behavior and collective excitations of the condensate. The validity of the Gross-Pitaevskii 
equation in describing the evolution of the condensate wave function has been verified by 
the quantitative comparison of theoretical predictions with experimental observations. 

The Gross-Pitaevskii equation (|) has been solved by means of a time-dependent varia- 
tional technique to obtain the low energy excitations of a trapped Bose gas for both positive 
and negative scattering lengths [3]. By taking a proper trial function of a fixed shape with 
some free time-dependent parameters, a variational technique based on Ritz's optimization 
procedure leads to a set of Newton-like ordinary differential equations of the second order 
for these parameters which characterize the solution to the NLSE. Avoiding expensive nu- 
merical simulations, this method affords analytical approximations which moreover provide 


a clear physical picture of the problem. 


For a dilute Bose gas trapped in a harmonic potential, with a Gaussian function taken 
as the natural trial wave function for low energy states, expressions for frequencies of the 
collective modes in а 3D anisotropic trap have been obtained [13]. In the large particle 
number limit, the expressions contain the spectrum derived by Stringari based on the 
Thomas-Fermi approximation, where the interatomic and trap interactions are dominant 
and the kinetic energy term is negligible. 

The finite-temperature excitation spectrum has also been measured in the experiments 
of Refs. | and [. The spectrum has been studied by means of the Popov version of 
the Hartree-Fock Bogoliubov approximation 1413) and by a two-fluid hydrodynamic ap- 
proximation ПОА. However, the finite temperature behavior of the spectrum observed in 
experiments has not been fully understood. In this paper, we extend the variational method 
at zero temperature to finite temperature where the thermal cloud exists. At a finite temper- 
ature, the thermal cloud provides an additional background field exerting on the condensate 
atoms, besides Ше mean-field interaction among them. Тће motion of the condensate can 
still be described by the nonlinear Schródinger equation with the mean-field modified to 
include the effect of the thermal cloud. The evolution equation for the condensate may be 
derived from a variational principle to minimize the action of the corresponding Lagrangian 
density 2. 


Let us consider а sample of Bosons confined in a d-dimensional harmonic potential 
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U(r) = Устава, (2) 
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where r = (z1,22,..., Фа). As shown in the recent experiment [fj], we may neglect the mode 


coupling between the condensate and the thermal cloud. In this simplication, we assume 


the static thermal cloud density distribution to be 


p»(r) = R(T) exp|-U (r)/kT], (3) 


where R(T) is the normalization factor. The condensate moves in the combined field of the 


external trapping potential and the mean-field potential due to both the condensate and the 


static thermal cloud. The Lagrangian density corresponding to the NLSE with the thermal 


cloud effect included can be written as 
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where the asterisk denotes a complex conjugate. The normalization condition for Ше macro- 


scopic wave function (г) is 
јачој = No, 


where Ло is the number of the condensate atoms. The thermal cloud density normalizes 


J arent) = NS 


where №, is the number of atoms in the thermal cloud. The last term in (Д) comes from the 


according to 


interaction between the condensate and the thermal cloud [18]. 


Following Ref. [13], we consider the Gaussian trial function 
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where, Бу setting the Gaussian distribution centered at the origin, we have ignored the irrel- 


evant motion of the center of the condensate. The variational parameters are the amplitude 
A, width w; and 6; relating to the curvature. All these parameters are real time-dependent 
functions, characterizing the macroscopic wave function of the condensate. Inserting the 
trial wave function into the Lagrangian density (И) and integrating over the coordinates, we 


obtain the Lagrangian 
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where we have used Ше equation of particle number conservation 


А0 = NL zm: 
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The last term in (1) is the contribution of the interaction between the condensate and Ше 
thermal cloud. 

By minimizing the action corresponding to the Lagrangian with respect to the parame- 
ters, we can derive the ordinary equations describing the evolution of the parameters. The 


widths of the condensate satisfy the following equation: 
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where, compared with the case of zero temperature, the thermal cloud contributes an extra 


term proportional to 
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In the derivation we have used the fact that the widths of the condensate are much smaller 
than those of the thermal cloud, i.e. mw?w;(t)? < kT. If we regard Ше widths (ил, wo, ...wa) 


as the coordinates of a fictitious particle, Eq. (B) is just Newton's equations of motion for 


the particle in the effective potential 
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where the w; ? term corresponds to the kinetic energy of the condensate, Ше f(T) term to 
the contribution from the interaction between the condensate and the thermal cloud, and 
the last a-dependent term comes from the interatomic interaction of the condensate atoms. 

After finding the widths, the curvature parameters 0;(Р) can be obtained through the 


following equation derived from the variation principle 
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Once the widths of the condensate are known, the rest of the parameters can be calculated, 
and the evolution of the Gaussian atomic cloud of the condensate can then be completely 
characterized. The whole problem of solving the NLSE reduces to that of solving the system 


of ordinary differential equations (В). 


As a special case, we first consider а two-dimension isotropic trap. For the ‘breath’ mode 


with a radial symmetry, Eq. (8) reduces to 
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where, for d — 2, 
f (T) = fo(T) = 81? N ah (w/2akTY’. 


After introducing the new dimensionless constant р = 1+2Май? / m? and new dimensionless 


variables 
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where ао = [h/(m4/1 — fo(T)w)]/?, Eq. (0) becomes 


v"(r) + v(T) = 


"OU (12) 


This equation, regarding v as the coordinate, describes the classical motion of a ficticious 
particle in the effective potential 

Veg —= v? T po, 
which admits an equilibrium point only when p > 0. In this case Eq. (12) has the analytical 
solution [19] 
E dde (à - z) zl . (13) 
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The oscillation frequency is 2074/1 — (Т). For a two-dimensional harmonic trap, №, ~ Т/? 


23], so the temperature dependency of the frequency is 


w(T) аа СТ, 


where C is a constant independent of T. At zero temperature the frequency reduces to 
2w, аз is expected. For a positive scattering length, the frequency decreases with increasing 
temperature. On the contrary, for a negative scattering length, the frequency increases with 


temperature. At a negative scattering length p becomes negative when No > та/(2јајћ“). 


This gives a critical value of the atom number beyond which по condensate would occur, 
nor is there any oscillatory solution. This critical value is independent of the temperature. 


In the case of a three-dimensional cylindrically symmetric harmonic trap 
L ee кой 
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the equations for the widths, in terms of the dimensionless constant Ро = 4/2/x No(T)a/ao 


and new dimensionless variables т = wt, и; = wj/ag with i = т, y, z and ag = 4/h/ mw, are 


du. 1 Po 
е yee 15 
dr? [t= (Ги v3 тм (a) 
и, 1 Po 
— + |1 — f(T = — + 16 
dr? | f3(T)] vy уз 7078 (16) 
Фи, Р 7 1 Po 
dr2 + [A n f(T) У; = y? + V2 Vy (17) 


Here, introducing the notation у = T/T., © = Aw and āo = J/Rhi/mo, we have denoted for 


The f(T) terms come from the interaction between the condensate and the thermal cloud. 


d=3 
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The v;* terms correspond to the kinetic energy contribution which is neglected in the 
Thomas-Fermi approximation but plays a role in stabilizing the condensate. As in the case 
of zero temperature, for a positive scattering length a > 0, the effective potential has a stable 
equilibrium point, while for a negative scattering length a < 0 there may be no equilibrium 
points, depending on No in comparison with a critical value N,. For an attractive interatomic 
interaction, Ne sets an upper limit of the atom number, i.e. a collapse occurs when № > №. 
When No « N,, a metastable equilibrium point exists as one of the two equilibrium points, 
supporting the condensate. In contrast to the two-dimensional case, for a three-dimensional 
trap, the critical atom number increases with temperature. This is consistent with what is 
known in the case of an isotropic trap [20]. 
Following Ref. [I3], we may study the collective excitation modes of the condensate Бу 


expanding и; of Eqs. (15)-(Г4) around the stable equilibrium point (vo, Мо, Voz), which satisfies 
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А linear analysis then leads to Ше following expressions for Ше low excitation frequencies 


о(т| = 2) = 2w/1 — p(T) -2Ра, (20) 
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where Ру = Ву (фи) and A*(T) = (X — f3(T))/(1 — fa(T)), and we have labeled the 
modes by Ше azimuthal angular quantum numbers m. Ву noticing that В(Г) = 0 at zero 
temperature, expressions (20) and (РТ) reproduce the results obtained in Ref. [I3]. 

For small Po, expressions (20) and (1) are close to the bare trap results. For large values 
of P) they correspond to the case of large interatomic interactions in comparison with the 
trap excitation energies. They are good for either the repulsive or attractive interatomic 
interactions, ie. for both the positive and negative scattering lengths. In the large atom 
number limit, where the kinetic energy term is negligible, expressions (РО) and (21) can be 


further simplified as 


w(|m| = 2) = wy 2(1 — fs(T)), (22) 


w(|m| = 0) = ш LAO) 4 + ЗА2(Т) + 4/16 +9М(7) — 1647) l (23) 


which, with (1— fa(t))w replaced by w and (А2 — fa(t))/(1— (ТУ) by A? at zero temperature, 
correspond to the Thomas-Fermi approximation. Near the critical temperature № is small, 
so the interatomic interaction of the condensate may not be dominant over the kinetic energy. 
For attractive interacting Bose gases of a negative scattering length, no condensate exists 
without the contribution of the kinetic energy being included. In these cases, instead of 
expressions (22) and (23), expressions (РО) and (21) should be used. 

In the MIT experiment [7], the axial trapping frequency и, = 16.93Hz and the radial 


frequency v, = 230Hz. Without considering the interatomic interaction, the critical tem- 
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perature kT, = hw(NA/1.202)'/3 [223] for N = 80 x 106 is 1.87иК, which is close to the 
experimental value T, = 1.7K. The difference may come from the neglected atomic mutual 


interaction [2324]. According to the local density approximation the chemical potential is 
given by 


pap, fg | Ш | = (7) | | | (24) 
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which, for very low temperatures, reduces to 
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Noticing that а = 2.75nm in the MIT experiment, we have the chemical potential и/Кв = 
348nk at No = 15 x 10° in agreement with the experimental value 380nK. This verifies the 
validity of the local density approximation. 

Using the experimental parameters we have calculated the oscillation frequencies. For 
the experiment on Rb vapor at JILA, the particle number is not too large, so expressions 
(Bd) and (B1) are used. The temperature dependent excitation spectrum from the theoretical 
calculation together with the experimental measurements is shown in Fig. 1. It is seen that 
at low temperatures the frequencies decrease with increasing temperature. A simple theory 
explains the experimental data at low temperature when the condensate occupies a signifi- 
cant proportion. In the MIT experiment for mode m = 0, the atom number is rather large, 
so we use expressions (22) and (23). The comparison of the theory with experiment is shown 
in Fig. 2. The theoretical results are in good agreement with experimental observation. 

In summary, by including a mean field acting on the condensate from the static thermal 
cloud in the Lagrangian density, we have studied the collective excitation modes of BEC 
at finite temperature by means of the time-dependent variational method, and obtained 
analytic expressions for the temperature-dependent frequency shift of low excitation modes. 
We have compared the theory with experimental observations. The theoretical results are 


in agreement with the MIT experiment very well. 
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FIGURES 


FIG. 1. Temperature-dependent excitation spectrum of the JILA experiment for collective 
modes m — 0 (triangle) and m — 2 (circle), compared with theoretical curves (dotted for m — 0 
and solid for m — 2). Frequencies are normalized with respect to the radial trap frequency, and 


temperature to the critical temperature for a harmonically confined ideal gas 


FIG. 2. Temperature-dependent collective excitation frequencies of the MIT experiment for 
mode т = 0 (circle), compared with the theoretical curve. In the experiment, temperature was 


varied by adjusting the rf frequency of a magnetic trap, so here Ам is a measure of temperature. 
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